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ABSTRACT 

The errors incurred in producing a thermal-inertia map are of three general 
types: measurement, analysis, and model simplification. To emphasize the geo- 

physical relevance of these errors, we express them in terms of uncertainty in 
thermal inertia and compare these uncertainties with the thermal-inertia values 
of geologic materials. Thus the applications and practical limitations of the 
technique are illustrated. 

All errors are calculated using the parameter values appropriate to a site 
at the Raft River, Id. Although these error values serve to illustrate the mag- 
nitudes that can be expected from the three general types of errors, extrapola- 
tion to other sites should be done using parameter values particular to the 
area . 


Measurement errors introduced by multispectral scanning systems commonly 
range from a noise-equivalent-temperature difference (NEAT) of 0.1K for air- 
craft systems to IK for satellite systems. The resulting uncertainties in 
thermal inertia range from 15 TIU (thermal- inertia unit)* to 150 TIU. 

Three surface temperature algorithms were evaluated: linear Fourier 

series, finite difference, and Laplace transform. In terms of resulting errors 
in thermal inertia, the Laplace-transform method is the most accurate ( 2 6 0 TIU), 
the forward finite-difference method is intermediate (300 TIU), and the linear 
Fourier series method the least accurate (460 TIU). However, the two more 
exact methods require more computer time, and both lack the ability of the 
Fourier-series algorithm to illustrate the physical significance of individual 
terms. By comparing the errors with the range of thermal inertias of geologic 
materials, it is possible to select the most cost/effective algorithm for a 
particular application. 

Model simplification errors result from three sources: transient effects, 

topography, and surface coating effects. Fo^ example, flux of 35 watts/m 2 
(equivalent to a water evaporation rate of 1.2 mm/day) would produce a thermal- 
inertia error of about 200 TIU. If no topographic corrections are made, a 10° 
southwestern slope causes an error in thermal inertia of about 350 TIU. A 
hematite surface coating one centimeter thick on a rock will produce an error 
of approximately 200 TIU, whereas a one-millimeter layer thickness will have a 
negligible effect (25 TIU). 

The total system errors in thermal inertia are placed in geologic context 
by noting the separation in thermal-inertia values for various geologic mate- 
rials. For example, the thermal-inertia separation between limestones and 
dolomites is typically 1200 TIU. The error analysis technique indicates that 
thermal-inertia discrimination should be possible between black shale and gabbro 
(separation of approximately 450 TIU); discrimination between these two mate- 
rials was not possible using Landsat spectral-reflectance data. 

Errors in thermal inertia can be translated into errors in bulk density 
and moisture content. Thermal-inertia mapping from aircraft (NEAT = 0.1K) has 
an accuracy of about 1 % in bulk density or equivalently, a sensitivity of 0.3? 
in water content. From satellite (NEAT = IK) these accuracies are 91 and 31 
respectively . 

A practical evaluation of the error analysis is demonstrated for aircraft 
data acquired at Raft River, Id. We cannot discriminate the tuff from the al- 
luvium, having a lower thermal inertia, and the lava, having a higher thermal 


TIU = 1 watt • secVm 2 /K 
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inertia, without a topographic correction, but we can discriminate the lava 
flows from the alluvium. 


1 . INTRODUCTION 

Thermal-inertia mapping is a recently developed technique that allows dis- 
crimination among geologic materials [1], [2], [3], [4], [5], [6]. The purpose 
of this paper is to discuss the types and magnitudes of errors encountered in 
applying this technique and to express these errors in a geologically meaning- 
ful way. The errors are of three major types: measurement errors, analysis 

errors, and model simplification and assumption errors. The major focus of 
this paper is to develop a procedure to select an algorithm which is appropri- 
ate to the geologic application and which minimizes computer time. 

In addition to presenting the mathematical basis for the three algorithms, 
we provide a numerical evaluation with respect to an exact solution. Estimates 
of the errors due to measurement and model simplification are combined with the 
analysis errors, and all errors are expressed as uncertainties in thermal 
inertia. These combined errors are compared to the thermal-inertia values 
determined for some selected geologic situations. As an illustration of the 
analysis technique, the Fourier-series algorithm is applied to the construction 
of a thermal-inertia map of an area in the Raft River, Id. 

2. MATHEMATICAL BASIS 

Thermal inertia is a physical property of geologic materials, and its value 
can be determined by measurements of the diurnal variation of surface tempera- 
ture and the albedo of the surface. A physical model for the surface tempera- 
ture variation is therefore required. 

We assume one-dimensional, periodic heating of a uniform half space of 
constant thermal properties. The ground temperature obeys the diffusion 
equation: 


3 2 v ( x , t ) _ 3 v(x,t) , . 

3x 2 3t 

where v(x,t) is the ground temperature at a distance x below the surface and 
a time t, and k is the thermal diffusivity. 

To solve this equation for the surface temperature, we require a form of 
the solution appropriate to the diffusion equation and the boundary condition, 
which assumes conservation of energy over the diurnal cycle. Before looking at 
possible forms of the solution, we need to investigate the boundary condition. 

If we assume a heat-balance condition exists at the surface, then the 
incoming heat fluxes — the incident solar radiation (I), the sky radiation (S), 
and the geothermal heat flux (Q) — must balance the outgoing fluxes — the con- 
vective flux (C), the evaporative flux (E), and the ground reradiation flux (R). 
The resulting expression of the heat flux at the surface is: 


=I+S+Q-C-E-R (2) 

x=0 

where K is the conductivity of the material. The convective and evaporative 
fluxes are included in equation 2 for completeness; however those terms will be 
handled only qualitatively in this paper. 

The parameter of interest, the surface temperature, is present in equation 
2 In the heat-conduction term (the left hand side of equation 2) and in the 
ground-reradiation term (R). Equation 2 is rewritten to display the surface 
temperature : 


3 v(x ,t ) 
3x 
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-K 


x = 0 


-Ea v(0,t) 4 + a(t) + Q 


( 3 ) 


3 vCxj.t) 
3x 


where e is the mean emissivity of the ground, o is the Stefan-Boltzmann con- 
stant, a(t) is a term combining the time-dependent solar flux and sky radiation 
flux and thus includes the surface albedo and emissivity, and Q Is the geo- 
thermal heat flux. This expression must be satisfied for any equational form 
of v(x,t) chosen for the solution to the diffusion equation. The nonlinear 
form of this boundary condition prevents a direct analytical solution of the 
diffusion equation. An exact numerical solution can be calculated [7]; however, 
the complexity, the difficulty of illustrating the physical significance of 
terms, and the computer costs all provide incentive for investigating less 
exact algorithms. 


One method of handling the nonlinear boundary condition is to linearize 
it and to express the solution to this condition in an exact form as an infi- 
nite Fourier series. The linearization is achieved by performing a Taylor-series 
expansion of the surface temperature around the sky temperature and discarding 
quadratic and higher order terms. This solution has been discussed in detail 
by Watson [l],.and the form of this solution is: 
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and T gk y is the sky temperature,. A is the ground albedo, C is an effective 

atmospheric-transmission factor, S p is the solar constant, Q is the geothermal 
heat flux, A^ and e n are the amplitude and phase of the n th harmonic of the 

insolation, and P is the thermal inertia. 


Another method of solving the diffusion equation is to use the exact, non- 
linear boundary condition and to apply a finite-difference iterative technique 
[9]. This solution assumes the form: 


( 2 (y _ V ) 

At • k m,n+l m,n 


( 5 ) 


where is the ground temperature at a time, nAt, and distance from the sur- 

face, mAx; Ax and At are small increments in distance and time; and k is the 
dif fusivity . 

The forward finite-difference method requires an initial solution and pro- 
pagates this solution forward in time. To insure continuity at the boundary, 
the heat-balance boundary condition must be satisfied. The boundary condition 
for the forward finite-difference solution is expressed as: 

V„ + l ' V l,n.l * IT f { V 0,n+l’ <n * 1> «} (6 > 

where K is the conductivity and f is the flux (eq. 2). 
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Numerical stability and convergence put constraints on the size of the 
time and distance increments. To insure these errors do not grow, the stabil- 
ity condition 

k • — < 0.5 (7) 

(At) 2 

must be met [ 9 ] . 

The third method of solving the diffusion equation is to use the nonlinear 
boundary condition and to apply Jaeger's method of Laplace transforms [7]. This 
method has the advantage that it directly yields the surface temperature with- 
out requiring an estimate of the temperatures at depths. However, the solution 
also must be iterated due to its nonlinear form: 

m 

P = .£= • ) ^ v <(> , n=l,2,...,m (8) 

n .rirr / > s *n-s+l * ’ ’ 

v s=l 

where F is the average flux into the ground in the ntti time interval; P is the 

thermal inertia; t is the period of the heating flux; v g is the average surface 

temperature in the sth time interval; and the <f>'s are a set of numerical coef- 
ficients determined solely by m, the total number of intervals in t. A detailed 
description of this method can be found in Jaeger [7] for the periodic heating 
of a half space, and the method was extended by Watson [8] to the heating of a 
layer over a half space. 


3. ANALYSIS ERRORS 

The previous section presents the mathematical basis for the three algo- 
rithms used to compute the diurnal surface temperatures (eq. 4, 5, 8). In this 
section, we numerically compare the errors by expressing the temperature dif- 
ferences as equivalent differences in thermal inertia. The difficulty in per- 
ceiving the physical significance of terms and the amount of computer time 
required for each algorithm are also discussed. 

The evaluation of the algorithms in terms of uncertainties in thermal 
inertia is performed in two steps. The first is to compare the three algorithms 
with each other by approximating the incident solar flux as a sinusoidal half 
wave. The second step is to compare the most exact algorithm — the Laplace 
transform — with an exact theoretical solution — a pure sinusoid flux. This 
establishes the relationships of the three methods to the exact solution. The 
resulting errors in thermal inertia for each of the algorithms compared to the 
exact solution are: 

Laplace-transform Algorithm 260 TIU 

Finite-difference Algorithm 300 TIU 

Fourier-series Algorithm 460 TIU 

Table 1 gives the list of parameter values used to make the above comparison. 

The relative ordering of the algorithms with respect to accuracy is as 
expected: the Laplace-transform algorithm, which is an implicit relationship 

between surface temperature and surface flux, is the most accurate; next the 
finite-difference solution; and lastly the Fourier-series algorithm, which 
linearizes the radiation-transfer terms. The magnitude of the errors can be 
made more geologically meaningful by noting that soils typically have thermal 
inertias in the range of 500 to 2000 TIU and that rocks typically have a range 
of thermal inertias of 2000 to 4000 TIU. (A more detailed discussion of the 
thermal inertias of geologic materials is presented in the discrimination of 
geologic materials section.) 

An important consideration in the cost/effectiveness of thermal-inertia 
mapping is the computer time required to solve the various algorithms. In the 
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above comparison, the finite-difference and the Laplace-transform algorithms 
were both iterated until a day/night surface-temperature difference of less 
than one degree was achieved between iterations. The Fourier-series algorithm 
requires the computation of a series summation but requires no iteration. The 
Laplace-transform algorithm requires twice as many iterations as does the 
finite-difference algorithm. The actual number of iterations depends upon the 
complexity of the boundary condition and the accuracy desired. Thus the rela- 
tive ranking of the algorithms in increasing order of required computer time 
is the Fourier series, the finite difference, and the Laplace transform. 

Another consideration in performing thermal-inertia mapping is the need to 
be able to see the relationship between various parameters and surface tempera- 
ture. From a knowledge of such relationships, the importance of various para- 
meters on the calculation of thermal inertia is seen. Only the Fourier-series 
solution (eq. 4) allows this awareness. For example, the form of the solution 
suggests a mean temperature level that is modulated by the amplitude and phase 
of the harmonic term. We can recognize that changes in the geothermal heat 
flux or changes in the emissivity will produce a shift in the diurnal tempera- 
ture curves; changes in the albedo will result in changes in both the mean level 
and the amplitude of the diurnal temperature. Thus a day/night temperature 
difference, which is a measure of this amplitude, does not contain, to first 
order, the effects of sky temperature, geothermal flux, or the zero-order term 
of the insolation. A recognition of these and similar relationships led to 
Watson [1] developing methods for both thermal inertia-mapping and geothermal 
heat-flux mapping. 

Table 2 summarizes the results of this section. 

1). TOTAL SYSTEM ERRORS 

Measurement errors introduced by the multispectral scanning system impose 
a fundamental instrumental limitation on thermal-inertia mapping. These errors 
are the result of a combination of factors, including the noise limit of the 
detector, amplifiers, tape recorder, and, for a satellite system, down-link 
telemetry. These errors are commonly expressed as noise-equivalent-temperature 
differences (NEAT). For an aircraft system, a typical value for a state-of-the 
art system is 0.1K. A IK NEAT is representative of a satellite system such as 
Nimbus III and IV; the HCMM (Heat Capacity Mapping Mission) satellite system 
may have a slightly lower NEAT. These values of uncertainty in temperature 
convert to approximately 15 TIU for the aircraft and 150 TIU for a satellite. 
These resultant errors are in addition to the errors introduced by the mapping 
technique and represent a fundamental limitation inherent to the technique. 

An error common to mapping data with both aircraft and satellite is that 
produced in the process of registering images acquired at different times. 
Because of changes in the vehicle path and scanner orientation (pitch, yaw, and 
roll), a point on the ground will not have the same image coordinates at two 
acquisition times. Daytime images acquired at the same time are automatically 
registered with each other and are geometrically identical. These "matched" 
images are called image-pairs. The registration is performed by selecting 
reference features on a master image-pair and subsequently identifying those 
features on the remaining "distorted" image-pairs and on the nightime thermal 
images. Reference features identified on both the master reflectance image and 
the master thermal Image must be identified on the "distorted" nighttime 
thermal images but may be located on either the reflectance or the thermal image 
of "distorted" daytime image-pairs. The choice depends upon whether the appear- 
ance of the feature is governed by topography, thermal inertia, or ref lectance . 
These reference features are then used as Input to a triangular interpolation 
algorithm . 

The errors due to misregistration are of two types. The first is a mis- 
registration of an image-pair to the master image-pair. This misregistration 
results in an association of incorrect albedos and/or incorrect temperature- 
difference computation — both of which result in an inaccurate calculation of 
thermal inertia. 


1151 


The second type of er - Is a misregistration of a calculated thermal 
inertia to a topographic m . This error is an inaccuracy in the assignment of 
a thermal inertia to a particular geographical point. We have not estimated 
the magnitude of these errors; however, we have qualitatively examined them. 

Misregistration of either type results from one of two conditions: the 

spatial frequency of the misregistration is higher than the spacing of the re- 
ference features or the image features are not discernible enough to select 
suitable reference points. The first of these conditions — the rapidly varying 
scene — is correctable by choosing more frequent reference points, and the re- 
sulting errors in thermal inertia can be minimized. The additional points do 
impose a significantly higher computer time. The second condition — low scene 
contrast — is common especially in nighttime imagery. In such situations the 
best estimate of the location of a reference feature on an image-pair is made 
by scaling the distance from other identifiable features. In this case, the 
resulting registration may be geometrically inaccurate but numerically equiva- 
lent due to the resulting low contrast image, which will not exhibit much vari- 
ation in temperature and/or reflectance values. Thus, these errors produced 
in computing or assigning thermal inertia are small. The edge effects that may 
be produced by identifying features by scaling can be reduced by spatial fil- 
tering . 

In addition to the measurement errors, another important source of error 
in thermal -inertia mapping results from simplifications introduced by a par- 
ticular model. These simplifications include the exclusion of transient fluxes, 
topography, and/or surface-coating effects. 

Transient effects are thos'e resulting from such conditions as cloudiness, 
windiness, or ground moisture. Effects of this type can be viewed as a flux 
perturbation on the periodic solution [11]. To illustrate this situation, we 
consider a heat-flux perturbation of 35 watts/m 2 . The magnitude of this per- 
turbation is equal to 2.51 of the solar constant, or an equivalent change in 
sky temperature of 10°K, or the effect of an evaporation rate of 1.4 x 10 - ® 
m/sec (1.2 mm/day ) . The transient-effect calculation was made by varying the 
duration of the perturbation and noting the resulting change in surface tem- 
perature at a later time. A graph expressing the resulting errors in thermal 
inertia is presented in figure 1. The maximum error introduced by this per- 
turbing flux is approximately 200 TIU. 

Correctable topographic effects require an assumption as to the scattering 
law for the sloping surface and the need to neglect reradiation between adja- 
cent surfaces. We shall assume that the surfaces are Lambert reflectors, and 
thus the local solar flux is proportional to the cosine of the angle between 
the surface normal and the solar radiation. In that case, an east or west 
slope to the surface causes a phase shift in the diurnal temperature curves, 
and a north or south slope to the surface produces a change in amplitude of the 
diurnal curves. A surface oriented in any other direction produces a combina- 
tion of phase shift and amplitude changes. Both changes result in erroneous 
day /night temperature differences. Similarly, slope and surface orientation 
affect the reflectivity and can cause erroneous determinations of albedo. Both 
temperature and albedo errors contribute to an incorrect determination of 
thermal inertia. The magnitude of this error (using the Raft River parameter 
values) for a 10° slope facing southwest is approximately 350 TIU with the max- 
imum error in thermal inertia occurring for a southwestern facing surface (fig- 
ure 2) . 

Surface coalings change both the surface reflectivity and the diurnal tem- 
perature; hence they affect the apparent thermal inertia of geologic materials. 
We examined the thermal effect of a hematite coating over a half space, having 
the average thermal inertia of igneous rocks. Figure 3 shows the implied 
change In equivalent thermal inertia for various thicknesses of this layer [8]. 
For example a hematite layer of 10 millimeters will produce a change in thermal 
inertia of approximately 200 TIU. If the layer th kness is less than 1 milli- 
meter, the coating Is transparent in terms of therr .-inertia mapping but will 
be observable as an albedo difference. 
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Of the three simplification effects described, only the topography cor- 
rection is at a stage of development to be incorporated into the thermal- 
inertia mapping. The local variations in the surface staining, vegetation 
cover, and surface roughness occur at ixscale which precludes detailed modeling 
because that would require detailed knowledge of these effects at every point 
in the test site. We consider these types of effects occurring at a local 
scale irreducible for any practical use of remote-sensing techniques. Similar- 
ly, topographic effects in localized areas due to departure from Lambert emis- 
sion, reradiation between adjacent surfaces, or shadowing effects are also 
probably in the same category. 

5. DISCRIMINATION OP GEOLOGIC MATERIALS 

In the previous sections, the types and magnitudes of errors encountered 
in performing thermal-inertia mapping were discussed. Now these errors are 
placed in geologic context and thereby provide criteria for evaluating the en- 
tire procedure. Thus the selection of modeling algorithm and data-collection 
procedure can be examined in the light of a particular geologic situation. 

As an illustration, assume that we wish to discriminate between limestone 
and dolomite from a satellite in a study site that is relatively flat and whose 
outcrops have iron-oxide stains of approximately one-millimeter thickness. 


Error Sources 


Error Magnitude 


satellite measurement 
topographic effects 
typical surface-coating effects 
typical transient effects 


150 TIU 
0 TIU 
25 TIU 
100 TIU 
275 TIU 


These measurement and model simplification errors are now combined with the 
analysis errors. The total system errors for each algorithm are: 


Algorithm 


Approximate Total Error 


Fourier series 
Finite difference 
Laplace transform 


750 TIU 
600 TIU 
550 TIU 


We now know that we cannot discriminate between geologic materials having 
a thermal-inertia difference of less than 750 TIU using the Fourier-series or 
of less than 550 TIU using the Laplace-transform algorithm. 

Figure M shows thermal-inertia histograms of various sedimentary rocks 
[12, 13]. From the histograms, we can determine that a typical thermal inertia 
for limestones is 2200 TIU and for dolomites is 3500 TIU; thus the separation 
is greater than the least accurate of the algorithms. Since any of the three 
analysis algorithms can be used, the selection of the algorithm will be deter- 
mined by the cost/effectiveness criteria and the Fourier series would be 
selected . 


In the hypothetical case, if we wish to discriminate between sandstones 
and dolomites, which typically have a thermal inertia separation of 550 TIU, 
we must select the Laplace-transform algorithm. 


Table 3 shows typical thermal-inertia values for a range of igneous rocks. 
The separation between a fine-grained felsic rock (rhyolite) and a coarse- 
grained felsic rock (granite) is only 200 TIU; the separation between rhyolite 
and a coarse-grained mafic rock (gabbro) is approximately M00 TIU. Discrimina- 
tion of a fine-grained felsic rock from a coarse-grained felsic rock can be 
achieved with the Laplace-transform technique, using many iterations and short- 
ening the time interval. Discrimination of a fine-grained felsic rock from a 
coarse-grained mafic rock can be achieved with either the finite-difference 
algorithm or the Laplace-transform algorithm, depending first upon the meas- 
urement and model errors and second upon the cost/effectiveness. 
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In the analysis of Landsat data, discrimination between black shale and 
gabbro was not possible using spectral-reflectance data [2]. The thermal 
inertia separation of these two materials is approximately 450 TIU (Table 4), 

[9, 13] and the magnitude of the separation indicates that discrimination by 
thermal inertia is practical even by satellite. Discrimination among basalt, 
andesite, argillite, and black chert was also reported to be not possible by 
Landsat. Table 4 shows that thermal-inertia discrimination should be possible 
between black shale and gabbro, between basalt/andesite and black chert, and 
between argillite and black chert. Discrimination among argillite, basalt, and 
andesite does not appear possible. 

Thermal inertia of soils is dependent primarily on bulk density and mois- 
ture content. Figure five shows that there is a rough linear relationship be- 
tween thermal inertia and density for dry materials. The addition of moisture 
to dry soil results in a rapid increase in thermal inertia. For example, dry 
sand has a thermal inertia of 590 TIU, and wet sand ( 8 % moisture) has a thermal 
inertia of 1020 TIU. Thus, for sand, the thermal-inertia mapping from aircraft 
(NEAT = 0.1K) has an accuracy of about 1 % in bulk density, or equivalently a 
sensitivity of 0-3? in water content. From satellite (NEAT = IK) these accu- 
racies are 9 % and 3 % respectively. 

6. RAFT RIVER ILLUSTRATION 

The method of analyzing errors in terms of uncertainties in thermal inertia 
was then applied to an analysis of a thermal-inertia map constructed of an area 
of the Raft River, Id. Reflectance and thermal data were acquired at five dif- 
ferent times in the diurnal cycle from July 24 to July 31, 1974. Figure six is 
the thermal image acquired near solar noon (July 31); figure seven is the ther- 
mal image acquired near solar midnight (July 31). These two images were used 
to form the temperature-difference image. Figure eight is the daytime reflec- 
tance image used to form the thermal-inertia map. 

Geometrically registered albedo and day/night temperature-difference 
images were formed according to the procedure outlined in the section on total 
system errors. A relative thermal-inertia image (fig -• 9) was then construct- 
ed from the images by forming the ratio of albedo to aay/night temperature dif- 
ference. The scale on the side of the image that relates the gray scale to 
absolute thermal-inertia values was determined from the Fourier-series algo- 
rithm using a sky temperature of 245°K and a cloud cover factor of 0.13. The 
sky temperature and cloud cover factor were determined using a least-squares 
model fit to all the repetitive thermal-scanner data for the site. In the for- 
mation of this image, no correction for topography was made. 

From the above information, the total measurement and analysis errors are 
compiled as follows : 

Errors Magnitude 

Fourier-series algorithm 460 TIU 

Aircraft measurement 15 TIU 

Three geologic materials have been discriminated on the thermal-inertia 
image [10]: lava flows, tuffs, and alluvium (figure 9). Histograms were con- 

structed of the thermal-inertia values of these materials from the map. The 
lava flows have a thermal inertia of 2010 + 250 TIU; the tuffs have a thermal 
inertia of 1470 + 210 TIU; and the alluvium has a thermal inertia of 1090 + 80 
TIU. From a topographic map, the alluvium was determined to be on a flat sur- 
face; the tuff area has a slope of 7-9° and an azimuth (measured counterclock - 
wise ) of 225 °; and the lava flows have a slope of 10° and an azimuth of 198°. 
Thus the error in thermal inertia due to topography is 0 TIU for alluvium, 

-95 TIU for tuffs, and +110 TIU for the lava flows. We therefore determined 
that we cannot discriminate the tuff from the alluvium and the lava without a 
topographic correction, but we can discriminate the lava flows from the 
alluvium. 
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7 . SUMMARY 

The utility of thermal-infrared surveys to discriminate geologic materials 
based on differences in their thermal properties has been recognized by pre- 
vious investigators. The error analysis technique presented here provides the 
means to evaluate the validity and reliability of such discrimination. This 
technique also provides criteria for selecting the most cost/effective algo- 
rithm for a particular geologic application. In addition, it can be applied to 
determine the fundamental accuracy required, and thus whether aircraft or sat- 
ellite (for example, HCMM) data would be appropriate. A parallel analysis 
technique could be developed for geothermal heat-flux mapping. Such a tech- 
nique could express the various errors in terms of the minimum detectable geo- 
thermal heat flux. 

Additional improvements in the mapping can be provided by incorporating a 
topographic correction into the model. One way of accomplishing this is to in- 
corporate digital topographic data in a look-up table [6]. Another method 
would be to assume a reflection law and compute the topographic correction 
using three reflectance images. Additional refinements would include the addi- 
tion of a diffuse illumination term in the insolation factor and a more detailed 
treatment of areal transient effects. In our treatment, the sky radiation and 
solar-transmission factor were assumed constant and were determined by a least- 
squares fit to the actual image data. A least-squares fitting including har- 
monic terms could be added to represent more accurately the actual areal 
transient effects. 

This study — even with its assumptions and simplifications on the model — 
should provide a useful comparison of the errors encountered in thermal-inertia 
mapping. Most importantly, all the uncertainties in the mapping technique are 
described in the context of discriminating geologic materials. 
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Figure 1. — Change in equivalent thermal 
inertia for a 35 watt/m 2 perturbing 
flux . 
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Figure 2. — Polar plot of change in 
equivalent thermal inertia for a 10° 

. slope for variable azimuth angles (e). 
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Figure 3* — Change in equivalent thermal 
inertia for a variable thickness hematite 
layer . 
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Figure 5- — Thermal inertia versus density for a variety of rock forming 
minerals, rocks, and soils. 
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Figure 7. — Thermal image acquired at 2300 hours of an area of the Raft River, 
Id. , July 31, 1974. 



Figure 8. — Reflectance image acquired at 1430 hours of an area of the Raft 
Fiver, Id. , July 31, 1974 • 


Figure 9- — Thermal-inertia map created from the 1430 and 2300 hours thermal 
image and a 1430 reflectance image of an area of the Raft River, Id. 
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TABLE 1. Parameter values used to compare the three algorithms. 



475 watts/m 2 

260°K 

1.0 


At = 8640 sec (0,1 day) 
Ax = 0.12 m 


TABLE 2. Relative ranking of algorithms with respect to cost/effectiveness. 


ALGORITHM 

Accuracy 

CRITERIA 

Computer Time Visualizing 

Physical Signi- 
ficance of Parameters 

Linear Fourier Series 

Least 

Least 

Yes 

Finite Difference 

Median 

Median 

No 

Laplace Transform 

Most 

Most 

No 


Table 3- Typical thermal-inertia values for a range of igneous rocks. 


ROCK TYPE 

THERMAL INERTIA 

(TIU) 

Rhyolite 

1950 


Granite 

2150 


Basalt 

2200 


Gabbro 

2350 


• 

4. Typical thermal-inertia vlues for 

some geologic materials that were 

studied using Landsat data [9, 13) 

• 


ROCK TYPE 

THERMAL INERTIA 

(TIU) 

Black Shale 

1900 


Gabbro 

2350 


Basalt /Andesite 

2200 


Argillite 

2250 


Black Chert 

3100 
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